Voltage fluctuations on a superconductor grain attached to a quantum wire 
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When a finite superconductor is in contact with a ID nor- 
mal conductor, superconducting phase fluctuations lead to 
power-law response of the normal subsystem. As a result, the 
charge fluctuations on the superconductor at zero tempera- 
ture have logarithmic correlator and a 1/w power spectrum 
(1//- noise). At higher temperatures 1/uj is pushed to the 
high frequency region, and 1/u 2 behavior prevails. 



Correlated character of transport in mesoscopic sys- 
tems strongly influences character of current and charge 
fluctuations there. Examples are non-Poissonian shot 
noise in quantum point contacts and NS junctions( ^] 
and references therein) and giant current noise in super- 
conducting quantum point contacts [^)-|J. Noise mea- 
surements may provide information about quantum cor- 
relations unattainable by any other means [^). 

In this paper we will show that equilibrium fluctua- 
tions in a mesoscopic NS system reveal the many-body 
reaction of a conducting system to a sudden external per- 
turbation, thus linking the noise to such phenomena as 
Fermi edge singularity (FES) || and Anderson orthog- 
onality catastrophe 0. Unlike the standard situation 
when the perturbation is due to creation of a scatterer, 
here the normal subsystem reacts to the change in the 
off- diagonal pairing potential in the superconductor, to 
which it is coupled via Andreev reflection processes on 
the NS boundary §,§. 

More specifically, consider the system shown in Fig.l. 
It consists of a mesoscopic superconducting grain S in 
contact with a ID normal wire N. We assume that the 
grain contains a large number of Cooper pairs, (n) 3> 1, 
and its capacitance C is large enough (so we are in the 
opposite limit to the Coulomb blockade regime [[l0| ) , but 
the contact region is smaller than the Josephson screen- 
ing length A,/, so that both the modulus and phase of 
the order parameter A = |A|e 1 ^ can be taken constant 
there. Temperature is low enough, T <C |A|, to neglect 
quasiparticle excitations in the superconductor. Finally, 
we assume that the ID conductor is scattering free, and 
its length L > ^o, where £o is the superconducting co- 
herence length. 

If isolated, the superconducting phase in the grain 
would fluctuate, the system being in the state with fixed 
particle number (see e.g. Ch.7). To the ID conduc- 
tor brought in contact with the grain, these phase fluctu- 



ations will translate into fluctuating Andreev boundary 
condition j8|. 

First let us present a heuristic picture of the phe- 
nomenon at zero temperature. The effect of a sud- 
den change in a boundary condition is a slow, power- 
law relaxation of the normal system to its new ground 
state, that is, Fermi edge singularity. It is conveniently 
described in terms of boundary conformal field theory, 
which allows a unified description of FES and Ander- 
son orthogonality catastrophe (llj. Both are character- 
ized by scaling dimension x of the boundary condition 
changing operator O, so that the correlation function 
(O^ (t)0(O)) ~ t~ 2x , and the dimension x itself is re- 
lated to the finite-size (0(1/ L)) correction to the ground 
state energy shift due to the action of O: 
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Hereafter we set fi = 1 and the velocity of the excitations 
(Fermi velocity in the non-interacting case) v = 1 unless 
stated otherwise; L is the characteristic size of the system 
(wire length in our case). This relation holds both for 
usual potential scattering and Andreev scattering || . 

The boundary condition changing operator in our case 
changes the phase of superconductor by $: 0$ = 
exp [3>jj^] = exp [j$n]. Here we used the relation between 
(j) and the Cooper pair number operator on the grain, 
" = 1 ~§4> > wn i cn holds if (n) ^> 1 JhJ . In agreement with 
general CFT argument, the Green's function of 0$ is 

G»(r) - (o* (r)0*(O)) = (e-^)^n^ = (t/t )- 2 * 

(2) 

with some cutoff parameter to. Though by itself G$ has 
no direct physical meaning, it is related to the correla- 
tor of charge fluctuations on the grain: Kq(t) / (4e 2 ) = 
([(5n(-r), (5n(0)] + ) /2, where Sh = n — (n) . Obviously, 



K q (t)-K q (0) 1 
Ae 2 2 



([6n(r),6n(Q)] + ) - (Sn 2 ) 



2 (K^U(°)J+) - (n )= hm p • 

Charge fluctuations are measured by e.g. monitoring the 
electrostatic potential of the grain. 

The dimension x is easily found if notice that the 
energy correction entering (Q) is simply the Josephson 
energy of a ballistic ID SNS junction of length L, with 
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phase difference $ ]12| , ^3[ at zero temperature: Ej{<&) = 
<J>7(4ttL), so that x = LEj{$)/tt = [$/(2tt)] 2 . Sub- 
stituting this and (0) into (^), we find for the charge 
fluctuations on the grain 



Kq(t) 
4e 2 
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(4) 



For the spectral density of charge fluctuations, Sq(u>) = 
2 J °° Kq(t) cos(wr)dr, such logarithmic behavior of cor- 
relator translates into Sq(uj) oc l/ui (1//- noise) in the 
corresponding frequency interval; it is often observed in 
spin glasses, magnetic polycrystals and other systems 
with wide distribution of relaxation times |l4|.|l5j . In our 
case the source of slow relaxation times is the normal 
system with its power-law approach to the new ground 
state. The idea of invoking infrared singularity to explain 
scaling that underlies current 1/ f- noise observed in nor- 
mal conductors was exploited in e.g. model of "quantum 
1/f-noise" due to Bremfitrahlung, but it encountered fa- 
tal difficulties (see |1| and references therein). In con- 
trast, we see no fundamental difficulty in the present 
mechanism, which will be shown to be consistent with 
fluctuation-dissipation theorem. 

Expression (|4|) is obviously only an intermediate 
asymptotics for 1 <C t/t <C exp[27r 2 (Sh 2 \]. Therefore 
we need a more consistent approach, taking into account 
effects of charging energy and finite temperature. We will 
see that the region of 1/f- fluctuations is pushed by ther- 
mal effects to high frequencies, and l/w 2 -type spectrum 
prevails except at very low temperature. 

We will describe the wire in terms of Tomonaga- 
Luttinger liquid (TLL), which allows us to treat simulta- 
neously non-interacting (ID Fermi-liquid) and interact- 
ing case of such S-TLL system |l6). Though as we will 
see interactions do not change qualitative picture of the 
phenomenon, they may lead to interesting side effects. 

Spin and charge degrees of freedom in TLL separate. 
Since spin is totally reflected by NS boundary both in 
the process of ideal Andreev reflection and ideal normal 
reflection, it is irrelevant and will be dropped. For the 
charge, the Lagrangian density is C = l/(2nK)[(dt9) 2 — 
(d x 9) 2 ], where K is the coupling constant, which is unity 
for non- interacting electrons. 

Physically, the boson field 9(x,t) corresponds to the 
phase of Charge Density Wave (CDW), while its dual 
<j>(x, t) to the superconducting phase of the Cooper pair 
||. For a perfect NS interface, the boundary value of 
the superconducting phase of the TLL should be equal 
to that of the superconductor, which in the absence of 
phase fluctuations gives a Neumann boundary condition 
for the field 9(x, t) in the wire |17| ], equivalent to a Dirich- 
let boundary condition on <f>: 0(0, t) = $ = const. 

The charge Q on the grain can be related to the bound- 
ary value of the "CDW" field 9. Indeed, the current in 
TLL can be written as 



J(x,t) 
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Therefore the electric charge of the superconductor at 
time t is given by 



Q(i) = Q(-oo) 



t fn 

J(0,t)dt = —6(0,t), (6) 



where we have defined Q and 9 so that Q = 9(0, t) = 
corresponds to the neutral superconductor grain. Simi- 
lar identification was used in studies of Coulomb Block- 
ade although our application is to quite different 
regime. 

Including the charging energy associated with Q, Eq = 

= ■j^p0(0, t) 2 , we arrive at the effective action of the 
problem: 



s = J J dtdx^ef - J dt^e(o,t) 2 



(7) 



Here 6 = 6/VttK and u = 2e 2 K/(irC) (u = 
2e 2 K/(nT%C) in conventional units). This is a free boson 
action with a mass term only on the boundary, and this 
effective theory is exactly solvable. The calculation of the 
charge fluctuation in the superconductor is reduced to 
the calculation of the correlation function of the bound- 
ary field 0(0, t). The physics of the problem is contained 
in a proper choice of the boundary condition at x = 
(on the NS boundary) @. 

The Matsubara Green's function of (for t = it) 
G{x,t) = - (T T e(x,r)e(0,0)>, is the Green's function of 
the Laplacian on the plane, with the boundary conditions 
d x Q(x = 0,r) = uQ{x = 0,t) and 9(x,t + /3) = Q(x,t), 
where (3 = l/T is the inverse temperature. As a conse- 
quence, the Green's function Q for the source at (#0,0) 
satisfies 



(d x 2 + d T 2 )G(x,T) =S(x-x )5{t) 
d x Q{x = 0,t) = uQ(x = 0,t). 



(8) 
(9) 



Because we are interested in the correlation function of 
the boundary, we take the limit xq — > 0. Integrating 
(^) over x from to xq + 0, we see that the boundary 
condition must be now modified: 



d x G(x = 0, r) = uG(x = 0, r) - 5{t). 



(10) 



For Fourier components we obtain an ordinary differen- 
tial equation: d 2 Q(x, ui n ) — uj n 2 Q with the boundary con- 
dition d x Q(x — 0, uj n ) — uQ(x — 0, Lu n ) — 1. Imposing the 
asymptotic condition lim^^oo Q(x, t) — > 0, we find 



G{x,u) n ) = 



-1 



u + \uj n \ 



(11) 



The Matsubara Green's function for the boundary field 
is Q[x = 0,w n ) = -l/(u + |w„|). 
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The correlation functions in the real time are obtained 
by an analytic continuation from the Matsubara Green's 
function. In a standard way, we obtain the retarded and 
advanced Green's functions, G R ( A '(w) = — l/(u^fiuj). 
Their mismatch at the real axis, which is pure imagi- 
nary, gives the spectral density of the system: p(u>) = 
— 2lmG R (uj) = 2uj/(u 2 + uj 2 ). The spectrum of electro- 
static potential fluctuations on the grain is given by: 



C 2 



^(Q(t)Q(0) + Q(Q)Q(t))e^ t dt 



2Ke 2 , ,1 , ,8u). 
-^pH-coth(-) 



2Ke 2 



irC 2 



— coth(^). (12) 



This result has the form expected from fluctuation- 
dissipation theorem (FDT) @. 

In the low temperature/high frequency limit (3u> 1, 
the spectrum is proportional to \oj\/(u 2 + uj 2 ). At fre- 
quency much higher than the "infrared cutoff' u, the 
spectrum is indeed l/u as expected from the heuristic 
argument (Q). On the other hand, at smaller uj fluctua- 
tions are suppressed, because finite capacitance obviously 
excludes large fluctuation of the charge in the grain. The 
l/u) spectrum, which is independent of the interaction 
parameter K. It automatically arises from the free boson 
field theory once the boson field itself is identified with a 
physical observable, because the correlation function of 
the boson field is logarithmic. We have thus shown an ex- 
ample of "quantum 1/f noise" naturally derived within 
the framework of equilibrium statistical mechanics and 
field theory. Although the l/io behavior is only expected 
at very low temperature and high frequency, we believe 
this presents certain interest both from the point of view 
of translation of infrared divergences into l//-spectrum, 
and of its destruction by thermal fluctuations. 

At higher temperature/lower frequency (3u> 3> 1, the 
spectrum is proportional to T/(u 2 + uj 2 ). The spectrum 
is (restoring the fundamental constants) AkTKe 2 / (irhui 2 ) 
for frequency larger than inverse capacitance uj u. In 
the opposite limit of small frequency to <C u, (but still 
u>/3 3> 1), the fluctuation is simply a white noise with the 
magnitude 2RnkT/(AK) where Rh = h/e 2 , again in a 
complete agreement with FDT. The transition from low- 
to high- temperature behavior occurs at [3uj ~ 1, namely 
kT ~ hf (see Fig.§). For T = lOmK, the boundary is 
at / ~ 200MHz. This would be within the capabilities 
of modern experiment; we expect the crossover between 
\/u) 2 and l/u) behavior could be directly observable. 

The charge fluctuation in the grain means the current 
flowing through the wire also fluctuates. The correla- 
tion function of the current can be obtained by a similar 
method. For example, the Matsubara Green's function of 
of the current at J(xq) is given by (2e 2 aj 2 /ir 2 )Gx {xo,uj n ), 
where 



Qx (xo,u n ) 



1 



2\U!„ 



1 - 



U + \UJ n 



is the Green's function of the Laplacian with the source at 
xq. The 1/uj fluctuation of the grain voltage translates to 
an oj-linear fluctuation of the current, although a similar 
behavior persists even for the current far from the grain. 

So far, we have considered a semi-infinite wire con- 
nected to the grain. In the remainder of this paper, we 
discuss how our results are affected by wire being finite. 

Instead of manufacturing a very long wire, it would 
be easier to connect the other end of wire to a normal 
reservoir. We assume that the junction to the reservoir is 
adiabatic (no backscattering) , and the phase information 
is immediately lost in the reservoir (perfect thermaliza- 
tion). Then the situation is not qualitatively different 
from the case of the semi-infinite wire, since the particles 
can leave to infinity (reservoir), and the usual reason- 
ing of Landauer formalism applies pjj ; the transport of 
the charge from/to the grain controlled by the geometric 
restriction (narrow ID channel) implies that our theory 
based on ID holds, even though the system is not purely 
ID. Thus we expect our results for the semi-infinite wire 
applies to the more realistic case of the finite wire con- 
nected to a normal reservoir, if the interaction effects in 
the wire can be neglected. 

When the interaction in the wire is not negligible, the 
problem is more subtle, since a quasiparticle excitation 
inside the wire is a collective motion of many electrons 
and cannot be simply transferred to the reservoir. While 
a rigorous theory of junction between a TLL and normal 
reservoir is not yet established, it has been studied by 
modelling the reservoir with a (semi-infinite) ID Fermi 
liquid. As a simplest model, we study a finite wire of 
length L, in which the coupling constant is K, attached to 
semi- infinite wire with coupling constant K' , K = aK 1 . 
To model the reservoir by a ID Fermi liquid, K' should be 
unity but we shall obtain the solution for general values 
of K' because it can be applied also to other interesting 
cases. It has been found that normal and Andreev-type 
reflections occur at the boundary if K ^ K ' [^2| . 

By a similar line of arguments as before, we obtain the 
Matsubara Green's function of the boundary field as 



5(0, 



[cosh (u> n L) + a sinh (|o; n |i)] 



(u + a\w n \) cosh (u> n L) + (\u n \ + au) sinh (ui n L) 



(14) 



(13) 



Since the temperature effect is given by the same FDT 
factor coth ((3uj /2)/2 as in the semi- infinite wire case, 
we concentrate on how the spectral density p(to) — 
-2lmG R (uj) is affected. 

For a = 1, the result reduces to the semi- infinite wire 
case, as it should be. For a « 1, the overall shape of 
the spectral density remains the same. However, a varia- 
tion with the period X/L is superimposed on the depen- 
dence for semi- infinite wire (Fig|j|). This is the resonant 
structure due to (normal and Andreev) reflections at the 
boundary between two wires. If the junction between the 



3 



wire and the reservoir is smooth and not abrupt, we ex- 
pect that the resonant structure is smeared out; then the 
observed spectrum would be quite close to the prediction 
for the semi-infinite wire. 

If a far from 1, the resonant structure dominates. In 
particular, the case K' — > (or K 1 — > oo) is equivalent to 
a finite wire with open end (or a finite wire attached to 
a very large superconductor reservoir), respectively. In 
these limits a — > oo, (a — > 0) the spectral density devel- 
ops into a set of discrete (5-functions with peak spacings 
of order 1/L, rather than a continuous function of u>. 
This reflects the finiteness of the system. In the L — > oo 
limit the (5-functions become dense and converge to the 
smooth function we have obtained for the semi-infinite 
wire. 

It is also instructive to take the opposite limit of short 
wire L — > 0. For the open end case (K' — > 0), the reso- 
nances are found at u ~ (2n — l)7r/ (2L) where n is a nat- 
ural number. For the wire connected to the superconduc- 
tor reservoir, in the short-wire limit, the resonances are 
found at to ~ tvk/L, n = 1,2,.... The resonance which 
would correspond to n = is actually affected by the fi- 
nite capacitance, and its location is given by w = 1/ \J~LC . 
This can be simply understood from the following ar- 
gument: from the conjugate relation between the phase 
and the charge, the effective Hamiltonian of the system 
is given by (see e.g. @) H = -^(^) 2 + JE^ ~ $ «) 2 > 
where $ is the superconducting phase of the grain, the 
first and second term represent the charging energy of the 
grain and the Josephson energy of the wire, respectively. 
This is a Hamiltonian of a harmonic oscillator with the 
eigenfrequency 1 / y/LC, exactly where the resonance is 
found. 

In conclusion, we have shown that quantum phase fluc- 
tuations in a superconducting grain in contact with a 
normal wire are translated into equilibrium charge fluc- 
tuations in the system, with specific spectrum demon- 
strating transition from 1/lu 2 to 1/uj behavior at very 
low temperature. The mechanism of fluctuations is re- 
action of the normal subsystem to sudden change of the 
off- diagonal pairing potential on the NS boundary (Fermi 
edge singularity) , leading in the low-temperature limit to 
quantum 1// noise. Similar effect can be expected to ex- 
ist also in systems with superconducting grain in contact 
with 2D conductors under appropriate conditions. 

While the experimental observation of this effect 
presents a certain challenge, we hope it to be possible 
in near future. Concerning the experiments, impurity 
scatterings and physical effects of probing will require 
further theoretical study, which is however beyond scope 
of the present paper. 
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FIG. 1. Superconducting grain S in contact with a ID con- 
ducting wire N. Tip P probes the electrostatic potential of 
the grain. 
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FIG. 2. Reduced spectral density of voltage noise in a 
grain, S(u) = S v (v)/[Ke 2 /nc 2 ] (Eq.(12)), as a function of 
Co = lo/u for different values of f s T/u. 
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FIG. 3. Resonant structure of spectral density for a = 0.1 
(dots), q = 1 (solid line), and a — 10. Dimensionless length 
of the wire L = 1; f = 1. Dashed straight line illustrates l/u> 
dependence. 
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